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The phylogeny of Japanese Cyperaceae was investigated using chloroplast ndhF sequence data from 37 
species in 27 genera. Our results show some congruence with the classification of Goetghebeur 
(1998). The 5.8S rDNA was analyzed for 37 species in 24 genera. Three base substitutions and 1-bp dele¬ 
tion were observed. A 3-bp indel (CAT) was found in many genera of Cyperaceae. Our results indicate 
that these mutations of 5.8S rDNA occurred in the ancestor of the family. The monophyly of subfam¬ 
ilies Mapanioideae, Sclerioideae, and Caricoideae is strongly supported. However, the monophyly of 
subfamily Cyperoideae is not supported because subfamily Caricoideae is nested within it. The 
monophyly of tribes Cariceae, Abildgaardieae, Eleocharideae, and Cypereae is strongly supported, but 
Schoeneae, Scirpeae, and Fuireneae are polyphyletic. Two groups were distinguished in the ndhF tree: 
(1) Dulichieae to Cariceae, and (2) Cypereae to Abildgaardieae. 
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The Cyperaceae are the third largest family 
of the Monocotyledons, comprising 104 genera and 
over 5000 species. They have a cosmopolitan dis¬ 
tribution, occurring in a variety of habitats ranging 
from the humid tropics to cold temperate regions 
(Goetghebeur 1998). In Japan, Ohwi (1944) recog¬ 
nized 20 genera and over 450 species. Koyama 
(1961), following the system of Ohwi (1944), rec¬ 
ognized four subfamilies and six tribes of Japanese 
Cyperaceae based mainly on floral morphology 
(Fig. 1). Ohwi (1944) and Koyama (1961) both 
hypothesized that the unisexual flowers of 
Cyperaceae evolved from bisexual flowers. The sex 
and morphology of the flower are therefore impor¬ 
tant characters for the subfamilial and tribal classi¬ 
fication of the family. 

The more recent systems of Bruhl (1995) and 


Goetghebeur (1998) were based on morphologi¬ 
cal, anatomical, embryological, phytochemical, and 
physiological analyses. Bruhl (1995) recognized 
two subfamilies and 12 tribes of Cyperaceae (Fig. 
1). His system was based on cladistic and phenetic 
analyses of DELTA (Description Language for 
Taxonomy) datasets (Dallwitz et al. 1993, Partridge 
et al. 1993). He classified four tribes with bisexual 
flowers in the subfamily Cyperoideae, and placed 
the other tribes into subfamily Caricoideae (Fig. 1). 
Goetghebeur (1998) recognized four subfamilies 
and 14 tribes, mainly based on floral and embryon¬ 
ic characters (Fig. 1). He classified nine tribes with 
bisexual flowers into two subfamilies Cyperoideae 
and Mapanioideae, and classified five tribes with 
unisexual flowers into two subfamilies Caricoideae 
and Sclerioideae. 
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Fig. 1 . Comparison of subfamilial and tribal classifications of Cyperaceae. 


Molecular data are providing a better under¬ 
standing of within- and between- tribal relation¬ 
ships in Cyperaceae (Simpson et al. 2003). Muasya 
et al. (1998) provided the first molecular phyloge¬ 
netic tree of the family using rbcL sequence data 
from 80 species in 40 genera. Muasya et al. (2000) 
revised the phylogeny of Cyperaceae using rbcL 
sequence and morphological data from 58 species in 
54 genera from Africa, Australia, North America, 
Europe, Southeast Asia, and South America. The 
molecular phylogenetic trees of Muasya et al. (1998, 
2000) do not support the subfamilial classifications 
of Bruhl (1995) and Goetghebeur (1998). Based 
on trnL-F region and rpsl6 gene (Simpson et al. 
2003), the phylogeny at the subfamilial level of 
Cyperaceae indicates that subfamily Mapanioideae 
is sister to the rest of the family. At present, only 


three molecular phylogenetic analyses have been 
conducted at the suprageneric level of Cyperaceae, 
none of which include the Japanese species (Muasya 
et al. 1998, 2000, Simpson et al. 2003). 

The rbcL data of Muasya et al. (1998, 2000) 
was insufficient for clarifying the relationships 
among the tribes of Cyperaceae. Other DNA data 
are necessary to provide a more detailed view of the 
relationships. Chloroplast ndhF has a higher rate of 
base substitution than rbcL, making it suitable for 
studies at the generic level and above (Olmstead & 
Sweere 1994, Clark et al. 1995, Bohs & Olmstead 
1997, Terry et al. 1997), therefore increasing the 
potential to better resolve suprageneric relation¬ 
ships in the family. 

Starr et al. (1999) reported a 3-bp insertion 
(CAT) near the 3' end of 5.8S ribosomal DNA 
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(rDNA) in Carex and Kobresia (position 155-157). 
Five types of 3-bp indels (CAC, CAT, CTC, CTT, 
and 3-bp deletion) were found in Eleocharis (Yano 
et al. 2004). They proposed that 3-bp indels are 
useful for delimiting clades within a genus or tribe. 
The mutations in the 3-bp indel, however, were 
investigated only in these three genera of 
Cyperaceae. 

The goal of this study is to clarify the subfa- 
milial and tribal relationships in Japanese 
Cyperaceae using chloroplast ndhF and nuclear 
5.8S rDNA sequence data, and to identify the 3-bp 
indels of 5.8 S rDNA in the Japanese members of the 
family. 

Materials and Methods 

The chloroplast ndhF gene from 37 species 
in 27 genera of Cyperaceae was sequenced and 
sequence data for three species were obtained from 
DDBJ/EMBL/GenBank (Table 1). Jimcus decipiens 
(Juncaceae) and Oryza sativa (Poaceae) were used 
as outgroups. Their ndhF sequences were deter¬ 
mined in this study or were obtained from 
DDBJ/EMBL/GenBank. The 5.8S rDNA region 
was analyzed for 37 species in 24 genera, using 
nine published sequences from Cyperaceae, three 
from Juncaceae and one from Oryza sativa 
(Poaceae) (Takaiwa et al. 1985, Roalson & Friar 
2000, Yano et al. 2004). 

DNA Extraction, PCR Amplification and DNA 
Sequencing 

Total DNA was extracted from O.lg fresh- 
frozen material or 0.03 g dried (silica gel dried or 
herbarium specimens) samples using the Nucleon 
Phytopure plant and fungal DNA extraction kit 
(Amersham Inc.). Total DNA preparations were 
purified by ethanol precipitation. The PCR reac¬ 
tion mixture contained 50 ng of genome DNA, 0.2 
pM of 5' and 3' end primers, 25 pL of Premix Taq 
(TaKaRa Bio Inc.), and 5% Dimethyl sulfoxide in a 
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total volume of 50 pL. Reactions were performed 
using a DNA thermal cycler GeneAmp PCR System 
2400 (Perkin Elmer Inc.). We designed four new 
primers: ndhF- 3F (5'-GGAACAGACATAT- 
CAATATG-3'), ndhF- 922F (5'-GAGCCACTT- 
TAGCTCTTGCT-3'), ndhF- 1163R (5'-TAAACC 
ACCCATAAGAACCA-3 ') and ndhF- 2016R (5'- 
TAACCACGATTATATGACCA-3'). Two primer 
sets were used for PCR amplification of the ndhF 
gene segments. One primer set was ndhF-3 F and 
ndhF-1163'R.; the other set was ndhF-9 22F and 
ndhF- 2016R. The reaction profile for ndhF ampli¬ 
fication involved: denaturation at 94° C for 5 min ; 40 
cycles at 92°C for 1 min, 48°C for 1 min and 72°C 
for 2 min; and a single final extension at 72° C for 7 
min. The nrDNA ITS fragments were amplified 
using two primers, ITSL (5-TCGTAACAA- 
GGTTTCCGTAGGTG-3'; Hsiao et al. 1994) and 
ITS4 (5 '-TCCTCCGCTTATTGATATGC-3'; White 
et al. 1990), following the protocol of Hsiao et al. 
(1994). 

The PCR products were electrophoresed on 
1.4% agarose gels to confirm a single band, and 
purified using the QIAquick PCR Purification Kit 
(QIAGEN Inc.). Direct DNA sequence analysis of 
the PCR products was performed using a DYEnamic 
ET Terminator Cycle Sequencing Kit (Amersham 
Biosciences Inc.) and analyzed using automated 
capillary electrophoresis ABI PRISM 310 (Applied 
Biosystems Inc.). The sequences obtained were 
submitted to the DDBJ/EMBL/GenBank databases. 

Sequence Analysis 

Data set - The species treated in this study 
are indicated in Table 1, along with their 
DDBJ/EMBL/GenBank accession numbers. For 
the analyses, we used two data sets. To estimate 
the phylogenetic relationships among the genera 
of Cyperaceae data set I (ndhF sequences) included 
39 species. The matrix of data set II (5.8S se¬ 
quences), including 41 species, was used to deter¬ 
mine the distribution of the unique 3-bp insertion in 
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Table 1. Taxa used in this study. Subfamilial and tribal classification of Cyperaceae according to Goetghebeur (1998). Vouchers are 
deposited in the herbarium of Okayama University of Science, Okayama (OKAY). Asterisks show literature citation for pre¬ 
viously published sequences. 


Taxon 

Cyperaceae 

Caricoideae 

Cariceae 

Carex albata Boott ex Franch. et Sav. 
Kobresia nepalensis (Nees) Kuk. 

Uncinia phleoides (Cav.) Pers. 

Cyperoideae 

Abildgaardieae 

Bulbostylis barbata (Rottb.) Kunth 

B. densa (Wall.) Hand. -Mazz. 

Fimbristylis velata R. Br. 

Cypereae 

Cyperus cyperoides (L.) O. Kuntze 

C. iria L. 

C odoratus L. 

Isolepis crassiuscula Hook. f. 

Kyllinga brevifolia Rottb. subsp. brevifolia 
Lipocctrpha microcephala (R. Br.) Kunth 
Pycreus sanguinolentus (Vahl) Nees 
Dulichieae 

Dulichium arimdinaceum (L.) Britt. 
Eleocharideae 

Eleocharis acicularis (L.) Roem. et Schult. 

f. longiseta (Svenson) T. Koyama 
E. congesta D. Don 

var.japonica (Miq.) T. Koyama 
E. interstincta (Vahl) Roem. et Schult. 

E. kuroguwai Ohwi 
E. mamillata Lindb. f. 

Fuireneae 

Actinoscirpus grossus (L. f.) 

Goetgh. et D. A. Simpson 
Bolboschoenus yagara (Ohwi) 

Y. C. YangetM. Zhan 
Fuirena umbellata Rottb. 

Schoenoplectus hotarui (Ohwi) Holub 
Schoeneae 

Cladium jamaicense Crantz 

subsp. chinense (Nees) T. Koyama 
Gahnia aspera (R. Br.) Spreng. 

G. tristis Nees 

Machaerina glomerata (Gaudich.) T. Koyama 
M. rubiginosa (Sol. ex G. Forst.) T. Koyama 
Rhynchospora chinensis Nees et Meyen 
R. rubra (Lour.) Makino 
Schoenus apogon Roem. et Schult. 

_ S. brevifolius R. Br._ 


Locality and voucher 


Accession number 
ndhF ITS 


Japan, Hokkaido Pref., Kyowa; Hoshino et al. 17400 
Nepal; Hoshino et al. 9660052 
Yen & Olmstead (2000)* 


AB261626 AB261661 
AB261627 AB261662 
AF163467 


Japan, Okayama Pref., Okayama ; Hoshino etal. 13779 
Japan, Okayama Pref., Soujya ;Hirahara etal. 19651 
Japan, Okayama Pref., Tamano; Hoshino et al. 17692 


AB261628 AB180718 
AB261629 AB261663 
AB261630 AB261664 


Japan, Okayama Pref., Yakake; Morimoto 17532 
Japan, Okayama Pref., Okayama; Yarn 18125 
Japan, Okayama Pref., Okayama; Katayama 19301 
Japan, Fukushima Pref., Shirakawa; 

Hirahara & Hoshino 19165 

Japan, Okayama Pref., Tamano; Hirahara & Yano 18135 

Japan, Okayama Pref., Sojya; Katayama 17518 

Japan, Okayama Pref., Okayama; Komagine & Masyo 17655 


AB261631 AB261665 
AB261632 AB261666 
AB261633 AB261667 
AB261634 AB261668 

AB261635 AB261669 
AB261636 AB261670 
AB261637 AB261671 


Yen & Olmstead (2000)* 

Japan, Nagano Pref., Karuizawa; Aral & Hirahara 18168 
Yano et al. (2004)* 

Roalson & Friar (2000)* 

Japan, Okayama Pref., Satosyo; Katayama 17643 
Yano et al. (2004)* 

Japan, Tokyo Pref., Hahajima; Katsuyama et al. 19915 

Japan, Aomori Pref., Misawa; Hoshino etal. 13020 

Japan, Okinawa Pref, Ishigaki; Hirahara & Yano 19320 
Japan, Hokkaido Pref., Tomakomai; Hirahara etal. 19269 

Japan, Okayama Pref., Ushimado; Izumi et al. 16881 


AF 163472 

AB261638 AB180693 

AB180697 

AF190611 
AB261639 AB 180704 
AB 180705 

AB261640 AB261672 

AB261641 AB261673 

AB261642 AB261674 
AB261643 AB206259 

AB261644 AB261675 


Japan, Tokyo Pref., Chichijima; Katsuyama et al. 19924 
Japan, Okinawa Pref., Iriomote; Hirahara 20361 
Japan, Tokyo Pref., Chichijima; Katsuyama et al. 19874 
Japan, Ehime Pref., Imabari; Hoshino et al. 19254 
Japan, Okayama Pref., Tamano; Hoshino et al. 17693 
Japan, Okayama Pref., Ushimado; Hoshino 18970 
Japan, Okayama Pref., Ushimado; Hoshino 18971 
Japan, Tokyo Pref, Chichijima; Katsuyama et al. 19875 


AB261645 AB261676 
AB261646 AB261677 
AB261647 AB261678 
AB261648 AB261679 
AB261649 AB261680 
AB261650 AB261681 
AB261651 AB261682 
AB261652 AB261683 
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Table 1. Continued. 


Taxon 

Locality and voucher 

Accession number 

ndhF ITS 

Cyperaceae 

Cyperoideae 

Scirpeae 

Eriophorum gracile K. Koch 

Japan, Hokkaido Pref., Tomakomai; Hoshino et al. 17382 

AB261653 

AB261684 

E. vaginatum L. 

Japan, Iwate Pref., Takizawa; Kameyama & Yano 18975 

AB261654 

AB261685 

Scirpus wichurae Boeck. 

Japan, Okayama Pref., Sojya; Katavama 17510 

AB261655 

AB261686 

Trichophorum alpinum (L.) Pers. 

Japan, Hokkaido Pref., Nakagawa; Sato 13260 

AB261656 

AB206270 

T. cespitosum (L.) Hartm. 

Japan, Iwate Pref, Ashiro; Katsuyama & Nakayama 20318 

AB261657 

AB261687 

Mapanioideae 

Hypolytreae 

Mapania paradoxa J. Raynal 

Givnish et al. (2002)* 

AY129256 


Sclerioideae 

Bisboeckelereae 

Diplacrum caricinum R. Br. 

Japan, Okayama Pref., Sojya; Katayama 19302 

AB261658 

AB261688 

Sclerieae 

Scleria parvula Steud. 

Japan, Okayama Pref., Sojya; Hoshino et al. 17555 

AB261659 

AB261689 

Outgroup 

Juncaceae 

Juncus decipiens (Buchenau) Nakai 

Japan, Okayama Pref., Okayama; Hirahara 17764 

AB261660 

AB261690 

Luzula capitata (Miq.) Miq. ex Kom. 

Japan, Okayama Pref., Okayama; Hirahara & Yano 17728 


AB261691 

L. elegans Guthnick 

Portugal; Draper 17729 


AB261692 

Poaceae 

Oryza sativa L. 

Hiratsuka et al. (1989)* 

X15901 



Takaiwa et al. (1985)* 


M16845 


the 5.8S gene (Starr et al. 1999) and to assess its 
phylogenetic utility. 

Alignment - DNA sequences were initially 
aligned using CLUSTAL W (Thompson et al. 1994) 
and then adjusted manually as necessary. 

Phylogenetic analysis 

Phylogenetic analyses were performed fol¬ 
lowing the modified methods of Tsubota et al. 
(2002, 2003) and Yano & Hoshino (2005). 
Phylogenetic trees were constructed using both the 
maximum-parsimony (MP) (Fitch 1971), and max¬ 
imum-likelihood (ML) (Felsenstein 1981) criteria. 
The trees obtained by the two methods were 
appraised by the log-likelihood measure with 
NucML in MOLPHY 2.3b3 (Adachi & Hasegawa 
1996). 

Parsimony analysis - In analysis with MEGA2 


(Kumar et al. 2001), the MP trees were searched by 
close-neighbor-interchange (CNI) with search level 
3, and the CNI search was repeated 10 times, start¬ 
ing with trees generated by 1000 replications with 
random addition sequences. In analysis by DNA- 
PARS in PHYLIP 3.573c (Felsenstein 1980-2001), 
the MP trees were searched for the best tree using a 
heuristic search. The parsimony method used the 
default setting. The sequences order was non ran¬ 
dom. Input sequences were sequential. In analysis 
by DNAPENNY in PHYLIP 3.573c, the MP trees 
were searched using branch-and-bound exact search¬ 
es. Again, the parsimony method used the default 
setting; the order was non random and input 
sequences were sequential. In analysis by 
PAUP*4.0bl0 (Swofford 2003), the MP trees were 
searched using the parsimony Ratchet search strat¬ 
egy (PAUPRat; Nixon 1999) in twenty 200-iteration 
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runs. 

Maximum likelihood analysis - ML trees were 
constructed by NucML in MOLPHY ver. 2.3b3, 
DNAML in PHYLIP 3.573c, fastDNAml ver. 1.2.2 
(Olsen et al. 1994). For the ML analyses, the 
HKY85 model (Hasegawa et al. 1985) was used as 
the estimate model. In analyses with NucML, the 
ML trees were searched with the local rearrange¬ 
ment method from the NJ tree by NucML and 
NJdist 1.2.5 (Adachi & Hasegawa 1996). In analy¬ 
ses by DNAML and fastDNAml, the ML trees were 
searched using global rearrangement options. The 
transition/transversion (TS/TV) parameter of 
DNAML and fastDNAml that best conforms to the 
data (4.644) was estimated based on calculations by 
MOLPHY. 

Tree comparison -Tree comparisons with the 
ML criteria by NucML was used to evaluate the 
topology of trees obtained from Consel O.lf 
(Shimodaira & Hasegawa 2001) and was also used 
to calculate the p- values of confidence for the bifur¬ 
cating candidate topologies using test procedures. 

Results 

The ndhF analysis 

The ndhF sequences of the 39 species, includ¬ 
ing two outgroups, ranged in length from 1893- 
2205 bp. Aligned sequences were 2248 bp, of which 
1716 sites were used for the phylogenetic analyses. 
Twenty-two gaps, ranging from one to twelve sites 
long, were introduced in the ndhF sequences, but 
gaps were treated as missing data. The data matrix 
of ndhF contained 761 variable sites, of which 484 
were potentially phylogenetically informative. 

For phylogenetic analyses based on ndhF 
sequences, a total of 24 topologies were obtained 
from the MP and ML methods. The 18 non-over¬ 
lapping topologies were obtained from the 24 
topologies by CONSEL O.lf. The best-supported 
ML tree was searched through the 18 non-overlap¬ 
ping topologies using a log-likelihood measure and 


^-values of AU, PP, and BP tests by CONSEL O.lf. 
The best-supported ML tree with the highest likeli¬ 
hood value is shown in Fig. 2. The log-likelihood 
value for the ML tree was -14229.79 + 353.5. AU, 
PP, and the BP tests values supported the ML tree 
more strongly than it supported the other trees. 

Our data provide strong support for Cyperaceae 
being monophyletic (100% Local bootstrap proba¬ 
bility: LBP, Fig. 2). Subfamily Mapanioideae is 
sister to the remaining species. Monophylies of 
subfamilies Sclerioideae and Caricoideae are sup¬ 
ported, respectively (99 and 100% LBP), but the 
subfamily Cyperoideae is not a monophyletic group, 
because the clade of Caricoideae is nested within 
Cyperoideae. Tribe Schoeneae is a polyphyletic 
group including three lineages. The first lineage is 
the genus Cladium and is sister to the rest of sub¬ 
family Cyperoideae and subfamily Caricoideae. 
The second lineage comprises three genera Gahnia, 
Schoenus and Machaerina, and each genus is mono- 
phyletic. The last lineage is the genus 
Rhynchospora, which is sister to the clade of the 
remaining Cyperoideae and Caricoideae. The tribes 
Cariceae, Scirpeae, Dulichieae, Abildgaardieae, 
Eleocharideae, Fuireneae, and Cypereae form a 
monophyletic group. This group comprises two 
large clades, Dulichieae to Cariceae and Cypereae 
to Abildgaardieae. Most tribes of this clade are 
monophyletic, but Fuireneae and Scirpeae are poly¬ 
phyletic. 

5.8S analysis 

The sequences of combined ITS 1 and 2 
ranged in length from 209 to 409 bp, and sequences 
of 5.8S rDNA were 163 to 166 bp, including 
Juncaceae and Poaceae. Of them, only sequences of 
5.8S rDNA were used, because the ITS regions 
were not strictly aligned. The data matrix of 5.8S 
(aligned length of 168 bp) contains 46 variable 
sites, of which 25 were potentially phylogenetical¬ 
ly informative and five were indels (Table 2). The 
variability in 5.8S rDNA sequences was too low 
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Table 2. The 25 informative sites and five indels in 5.8S rDNA. Shaded boxes show the noticeable mutation. 

000000000000111111111111111111 

000011123369002233333333555666 

123512870470138901234568567478 


Cyperaceae 
tribe Cypereae 

Kyllinga brevifolia subsp. brevifolia 
Cyperus odoratus 
Pycreus sanguinolentus 
Cyperus iria 
Cyperus cyperoides 
Lipocarpha microcephala 
Isolepis crassiuscula 
tribe Fuireneae 

Schoenoplectus hotarui 
Actinoscirpus grossus 
Fuirena umbellata 
Bolboschoenus yagara 
tribe Eleocharideae 
Eleocharis mamillata 
Eleocharis congesta vix.japonlca 
Eleocharis interstincta 
Eleocharis kuroguwai 
Eleocharis acicularis f. longiseta 
tribe Abildgaardieae 
Bulbostylis densa 
Bulbostylis barbala 
Fimbristylis velata 
tribe Scirpeae 

Trichophorum alpinum 
Trichophorum cespitosum 
Scirpus wichurae 
Eriophorum vaginatum 
Eriophorum gracile 
tribe Cariceae 

Kobresia nepalensis 
Carex albata 
tribe Schoeneae 

Rhynchospora chinensis 
Rhynchospora rubra 
Gahnia aspera 
Gahnia tristis 
Schoenus apogon 
Schoenus brevifolius 
Machaerina rubiginosa 
Machaerina glomerata 
Cladium jamaicense subsp. chinense 
tribe Sclerieae 
Scleria parvula 
tribe Bisboeckelereae 
Diplacrum caricinum 
Juncaceae 

Juncus decipiens 
Luzula elegans 
Luzula capitata 
Poaceae 

Oryza sativa 


AGTTTTCCGGGTGGA-AGGG-TGG-T A A 

ATTCTCCCCCGTGCG-ACCC-TCC-T A A 

ATTTTCCCCCGTGCG-ACCC-TCC-T A A 

AATTTCCCCCGTGCG-ATCC-TCC-T A A 

AATTTCCCCCGTGCG-ACCC-TCC-T A A 

AATTTCCCCCGTGCG-ACCC-TCC-T A A 

ATTTTCCCCCGTGCG-ATCC-TCC-T A A 

AATCTCCCCCGTGCG-ATCC-TCC-T A A 

ATTTTCCCCCGTGAG-ACCC-TCC-T A A 

TGTCTCCTCCATGCG-ACCA-TCCCATTAA 

AATTTCCCCGGTGCG-ACCC-TGCCATTAA 

GAGT CCCCCCATGCG-ACCA-TCC-T A A 

GAGTCCCCCCATGCG-ACCA-TCCCTTTAA 

GAGCCCCCCCATGCG-ACCA-TCCCTCTAA 

GAGTCCCCCCATGCG-ACCA-TCCCACTAA 

GAGTCCCCCCATGCG-ACCA-TCCCATTAA 

GAGTCCCCCCGTGCG-ACCA-TCCCATAAA 

GAGTCCCCCCGTGCG-ACCA-TCCCATAAA 

GAGTCCCCCCGTGCG-ATCA-TCTCATAAA 

TATCTCCCCCGTGCG-ACCC-TCCCATTAA 

TATTTCCCCCGTGCG-ACCC-TCCCATTAA 

AGTTTCCCCCGTGCG-ACCC-TCCCATTAA 

AATTTCCCCGGTGCG-ACCC-TCCGATTAA 

AATCTCCCCCGTGCG-ACCC-TCCCATTAA 

AGTTTCCCCCGTGCG-ACCC-TCCCATTAA 

AATTTCCCCCGTGCG-ACCC-TCCCATTAA 

AAGTTCCCCCGTGCG-ACCC-TCCCATTAA 
AAGTTCCTCCGCGCG-ACCC-TCCCATTAA 
AT GTTCCCCCGCGTG-ACCA-TCCCATTAA 
ATGTTCCCCCGCGCG-ACCA-TCCCATTAA 
AT GTTCCCCCGCGCG-ATCC-GCCCATTAA 

ACGTTCCCCCGTGCG-ATCC-GCC-T A A 

TTGTTCCCCCGCGCG-ACCC-GCC-T A A 

TCGTTCCCCCGCGCG-ACCT-GCC-T T A 

AT GTTCCCCCGCGCG-ACCC-GCCCATTGA 

AACTTCCCCCGCGTG-TCCT-GCCCATAAA 

GAGTTCCCCCGTGCG-GCCA-ACCCATTAA 

GGATTCCCCCGTGTCCATTT-GGT-CTT 

GAAATTTTTTGTATCCACTTTGGT-CTT 

GAACTCTCTTGTATCCTTTTTGGC-CTT 

CACCTCCCCCCTGCCCATCC-GGC-CGC 
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Fig. 2. The highest likelihood tree for the aligned 1716 bp of 39 ndhF gene sequences (HKY 85 model; alp — 4.644; In L = -14229.79 
± 353.5 by NucML). The subfamilial and tribal classification by Goetghebeur (1998) is shown on the right. The root is arbitrarily 
placed on the branch leading to Mapania paradoxa. Local bootstrap probabilities (LBP; %) are shown near branches. The taxa 

were not investigated of 5.8S rDNA. Three large capitals show a 3-bp indels of 5.8S rDNA 3-bp deletion). Bars across 
branches are supported by the noticeable mutations of 5.8S rDNA <the sequence number>. 


to resolve the phylogenetic relationship at the supra- 
generic level. Eight sites were noticeable mutations 
arranged on the tree (Table 2, Fig. 2). Three base 
substitutions and 1-bp deletion (positions 128, 129, 
136, and 168) were observed in Cyperaceae. Of 
the substitutions, for the single base substitution at 
position 128, the genus Kyllinga shows adenine 
(A), but the other genera show guanine (G). A sin¬ 


gle base substitution (position 11) is observed in 
tribes Eleocharideae and Abildgaardieae. A sin¬ 
gle base substitution (position 135) is shared in 
tribes Schoeneae (excluding Machaerina, Schoenus, 
and Cladium), Cariceae, Scirpeae, Dulichieae, 
Eleocharideae, Abildgaardieae, Fuireneae, and 
Cypereae. In the same position, another single sub¬ 
stitution is observed in Diplacrum of tribe 
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Bisboeckelereae. A 3-bp insertion (CAT; position 
155-157) is newly observed in the 12 genera of 
tribes other than tribe Cariceae. Eight genera of 
Cyperaceae do not have a 3-bp insertion. In addition, 
two types of 3-bp indels (CAT and 3-bp deletion) 
are found in Schoenus. Five types of 3-bp indels 
(CAC, CAT, CTC, CTT, and 3-bp deletion) were 
ascertained in Eleocharis (Table 2, Yano et al. 
2004). 

Discussion 

The ndhF phylogeny 

This study represents the first comprehensive 
phylogenetic analysis of Japanese Cyperaceae using 
ndhF. The ndhF analysis of suprageneric relation¬ 
ships in the family has shown better resolution than 
the rbcL analyses of Muasya et al. (1998, 2000). 
Our analyses support the tribal classification of 
Goetghebeur (1998), but the monophyly of sub¬ 
family Cyperoideae is not supported. Our data also 
show no congruence of the subfamilial classifica¬ 
tions of Koyama (1961) and Bruhl (1995). 

Mutations of 5.8S rDNA 

Three base substitutions and 1-bp deletion 
(positions 128, 129, 136, and 168) support the 
monophyly of Cyperaceae (Table 2, Fig. 2). Starr et 
al. (1999) reported a 3-bp insertion (CAT) near the 
3' end of the 5.8S region in tribe Cariceae (position 
155-157). Yano et al. (2004) found five types of 3- 
bp indels (CAC, CAT, CTC, CTT, and 3-bp dele¬ 
tion) in Eleocharis. In our analysis, a 3-bp indel 
(CAT) is found in many genera of Cyperaceae, and 
3-bp deletion is found in 10 genera (Table 2). Our 
results indicate that a 3-bp insertion (CAT) occurred 
in the ancestor of Cyperaceae, and a 3-bp deletion 
occurred in four different evolutionary lineages 
(Fig. 2). Starr et al. (1999) and Yano et al. (2004) 
proposed that 3-bp indels are useful to delimit a 
clade within a genus or tribe. In our analyses, a 3-bp 
deletion supports the Kyllinga to Actinoscirpus and 
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Machaerina clades. In addition, Schoenus is divid¬ 
ed into two groups by 3-bp indels (CAT or 3-bp 
deletion). These results suggest that 3-bp indels are 
useful in obtaining a better understanding of intra- 
and inter- generic relationships in Cyperaceae. 

Subfamily Mapanioideae 

Goetghebeur (1998) recognized four subfam¬ 
ilies in Cyperaceae: Mapanioideae, Sclerioideae, 
Caricoideae, and Cyperoideae (Fig. 1). Subfamily 
Mapanioideae was also recognized by Koyama 
(1961). It is generally recognized that the 
Mapanioideae are the most primitive subfamily of 
Cyperaceae, based on floral morphology (Holttum 
1948, Kern 1974). In the ndhF tree, subfamily 
Mapanioideae is placed as a sister to the rest of 
the Cyperaceae (Fig. 2). A similar position for sub¬ 
family Mapanioideae was found by Muasya et al. 
(1998, 2000) and Simpson et al. (2003). Our find¬ 
ings support the belief that subfamily Mapanioideae 
is the most primitive subfamily in Cyperaceae. 

Subfamily Sclerioideae 

Subfamilies Sclerioideae and Caricoideae have 
spicate spikelets composed of unisexual flowers. 
However, subfamily Caricoideae differs from sub¬ 
family Sclerioideae by having perigynia in the pis¬ 
tillate flowers (Koyama 1961, Goetghebeur 1998). 
Goetghebeur (1998) recognized four tribes 
Bisboeckelereae, Cryptangieae, Sclerieae, and 
Trilepideae in subfamily Sclerioideae. In the ndhF 
tree, tribes Bisboeckelereae and Sclerieae form a 
clade (99% LBP). Muasya et al. (1998, 2000) also 
reported tribe Bisboeckelereae to form a clade with 
tribe Sclerieae using rbcL analyses (93 or 99% BP: 
bootstrap probabilities). 

Subfamilies Cyperoideae and Caricoideae 

Goetghebeur (1998) merged seven tribes with 
spicate spikelets composed of bisexual flowers into 
subfamily Cyperoideae, and classified tribe 
Cariceae, which have unisexual flowers and a perig- 
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ynia, into subfamily Caricoideae. Several morpho¬ 
logical analyses (Bentham 1883, Clarke 1908, 
Mattfeld 1936, Koyama 1961, Bruhl 1995) sug¬ 
gested that subfamily Caricoideae are closely relat¬ 
ed to tribe Sclerieae, which were treated as a mem¬ 
ber of subfamily Sclerioideae in Goetghebeur 
(1998). In the rbcL trees of Muasya et al. (1998, 
2000), subfamily Caricoideae is not related to tribe 
Sclerieae, but is nested within the clade of sub¬ 
family Cyperoideae with low bootstrap support. 
Our ndhF tree strongly supports the monophyly of 
subfamily Caricoideae and its nesting within sub¬ 
family Cyperoideae. It also indicates that subfami¬ 
ly Caricoideae is related to tribe Scirpeae of sub¬ 
family Cyperoideae with strong support (100% 
LBP). Taxonomic revision is therefore needed for 
subfamilies Cyperoideae and Caricoideae. 

Based on our findings, we propose that subfamily 
Caricoideae be treated as a member of subfamily 
Cyperoideae. 

Phylogenetic relationship among tribes in the clade 
Cyperoideae 

Goetghebeur (1998) classified seven tribes, 
Abildgaardieae, Cypereae, Dulichieae, 
Eleocharideae, Fuireneae, and Scirpeae, in sub¬ 
family Cyperoideae (Fig. 1). All tribes except 
Schoeneae form a monophyletic group with the 
Cariceae in subfamily Caricoideae in the ndhF 
tree, with two clades distinguished within this group: 
(1) Dulichieae to Cariceae, and (2) Cypereae to 
Abildgaardieae (Fig. 2). Tribes in the first clade 
are mainly distributed in the northern hemisphere; 
tribes in the second clade have a cosmopolitan dis¬ 
tribution (Goetghebeur 1998). Our data show that 
these two clades diverged early in the evolution of 
the Cypereae to Cariceae clade. 

The second clade comprises of two subclades, 
Eleocharideae to Abildgaardieae and Cypereae to 
Fuireneae. Svenson (1939) proposed that Eleocharis 
of tribe Eleocharideae is closely related to the 
Fimbristylis of tribe Abildgaardieae based on mor- 
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phological characters. Recent morphological stud¬ 
ies (Tucker 1987, Bruhl 1995, Goetghebeur 1998), 
however, proposed that the Eleocharideae are close¬ 
ly related to Fuireneae. Muasya et al. (1998, 2000), 
based on by rbcL analyses (< 50 or 54% BP), 
reported that tribes Eleocharideae and Fuireneae 
formed a clade. Roalson & Friar (2000) and Yano et 
al. (2004), however, reported tribes Eleocharideae 
and Abildgaardieae form a clade (52 or 86% BP) 
based on ITS analyses. In the ndhF tree, tribes 
Eleocharideae and Abildgaardieae form a clade 
(93% LBP). In the 5.8S rDNA, a single base sub¬ 
stitution is shared in these two tribes (position 11). 
Our data support a shared more recent common 
ancestry of tribes Eleocharideae and Abildgaar¬ 
dieae, but not Fuireneae. 

In the tribe Schoeneae, three clades were dis¬ 
tinguished by ndhF analysis (Fig. 2). In addition, 
tribe Schoeneae was divided into two groups by a 
single base substitution in the 5.8S rDNA (posi¬ 
tion 135). Our data indicate that tribe Schoeneae is 
polyphyletic. Taxonomic revision is needed for 
tribe Schoeneae. 
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